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ABSTRACT 

A variety of modifications to the modeled dissipation rate transport equation that have 
been proposed during the past two decades to account for rotational strains are examined. 
The models are subjected to two crucial test cases: the decay of isotropic turbulence in a 
rotating frame and homogeneous shear flow in a rotating frame. It is demonstrated that 
these modifications do not yield substantially improved predictions for these two test cases 
and in many instances give rise to unphysical behavior. An alternative proposal, based on 
the use of the tensor dissipation rate, is made for the development of improved models. 
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1. Introduction 


Many turbulent flows of practical importance involve combinations of rotational and 
irrotational mean strains. The rotational strains can arise from the effects of curvature, a 
system rotation, or swirl, as well as from shear. Much of the motivation for formulating 
second-order closure models has its origins in the desire to account for the distinct physical 
effects of rotational and irrotational strains, in a framework where directional information 
is incorporated. Reynolds stress models of the eddy viscosity type, since they only depend 
on the mean velocity gradients through their symmetric part (i.e., the mean rate of strain 
tensor), are oblivious to the presence of rotational strains. Consequently, these models 
cannot distinguish between the considerably different turbulence physics exhibited in plane 
shear, plane strain, and rotating plane shear flows - a deficiency that is now well-known 
[1, 2]. However, the same deficiency still remains with the modeled form of the dissipation 
rate transport equation that has been commonly used in the turbulence modeling literature 
starting with the work of Launder and co-workers [3]. This frequently used model has 
no explicit dependence on rotational strains; rotational strains can only have an indirect 
effect through the changes that they induce in the Reynolds stress tensor. Such a limited 
dependence is known to be defective, making it impossible to properly describe rotating 
isotropic turbulence and swirling jets, among other turbulent flows [4, 5]. Consequently, 
several independent researchers over the past two decades have tried to develop rotational 
modifications to the modeled dissipation rate transport equation to resolve this problem. 
These alterations will be the subject of the present study. 

In this paper, a variety of previously proposed modifications to the modeled dissipation 
rate transport equation will be considered: 

(1) the Pope [5] model, 

(2) the Hanjalic and Launder [6] model, 

(3) the Bardina model [7], and 

(4) the Raj [8] model. 

These models represent a good cross-section of the rotational modifications to the modeled 
dissipation rate transport equation that have been proposed during the last fifteen years. 
Two flows from homogeneous turbulence will be used to test the models: (a) the decay of 
isotropic turbulence in a rotating frame, and (b) homogeneous shear flow in a rotating frame. 
These two test cases are chosen since they represent homogeneous flows whose structure is 
significantly altered by rotation, without the added complication of walls or other sources 
of mean turbulent diffusion. In this fashion, the question of the interaction of rotational 
and irrotational strains can be studied in isolation, independent of the complicating features 
introduced by solid boundaries and mean turbulent diffusion. The performance of each of the 
models will be documented in detail and specific proposals will be made for the development 
of improved models. A new model based on the tensor dissipation - formulated in order to 
include more directional information - will be proposed and tested. 
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2. The Models and Test Cases to be Considered 


We will restrict our attention to the incompressible flow of a viscous fluid with constant 
properties. The velocity field v and pressure P are decomposed into mean and fluctuating 
parts, respectively, as follows: 


v = V + u, P = P + p. 


( 1 ) 


The exact transport equation for the turbulent dissipation rate in homogeneous flows takes 
the form [2] 

£ = P'-$' (2) 

where 

_ duj dukdvj _ „ duj du } dv t _ du t du t du fc ~ 

e dxj dxj dxk dx{ dxk dxk dxk dxj dxj 


= 2v 


2 d 2 Ui d 2 Ui 


(4) 


dxjdxk dxjdxk 

ar e, respectively, the production and destruction of dissipation terms given that e = 
2i/dv,i / dxj dui j dxj is the scalar dissipation rate and v is the kinematic viscosity of the fluid. 
Here, since the turbulence is homogeneous, the overbar denotes a spatial mean, and the 
superposed dot denotes an ordinary time derivative. 

The modeled version of the transport equation for e that has been used widely in the 
turbulence modeling community - which is due to Launder and co-workers [3] - takes the 
form 

(5) 


e = C tl ^V-Cj- 


where V = —uiu]dv t /dxj is the turbulence production, and C e \ and C e 2 are dimensionless 
constants that typically assume the values of 1.44 and 1.92, respectively. Equation (5) 
is obtained from (2) by making the assumption that the production (or destruction) of 
dissipation is proportional to the production (or destruction) of turbulent kinetic energy, 
i.e., 

V e oc V/t 0 , oc e/r 0 (6) 

where To is the turbulent time scale that is needed for dimensional consistency. The added 
assumption that the turbulent time scale r 0 is given by 


K 


TO = 


(7) 


then yields the modeled equation (5). 

Due to the symmetry of the Reynolds stress tensor = u t Uj, it follows that 


V = -T^Sij 


( 8 ) 


where 


1 / dfJj dvj 

2 \dxj dxi j 


(9) 
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is the mean rate of strain tensor. Hence, (5) has no explicit dependence on the mean rotation 
tensor 

1 

Ui-i — — 

tJ 2 

It therefore follows that the commonly used modeled dissipation rate equation can only 
be affected indirectly by rotational strains through the changes that they induce in the 
Reynolds stress tensor. As we will show later, this makes it difficult to properly describe 
certain rotating turbulent flows. 

The model suggested by Pope [5] is based on the addition of a mean vortex stretching 
term and takes the form 

£ = C'el’^'P — + C e3 —^SijlJjkU}ki (11) 

in homogeneous turbulence where C e i = 1.44, C c 2 = 1.92, and C e3 — 0.79. The standard 
e-equation is obtained in the limit as C e3 — * 0. In the modification proposed by Hanjalic 
and Launder [6] to “sensitize the dissipation equation to irrotational strains,” the e-equation 
takes the following form in homogeneous turbulence: 

* = C ' l ~K V ~ Cc2 K ~ C * Ke ™i e ™ kl dx . Q Xl ( 12 ) 

where is the permutation tensor and C c i = 1.44, C e 2 = 1.92, and C c 3 = 0.27. Since, 

+ (13) 

it is clear that this modification introduces rotational strains; in the limit as C e3 — > 0 the 
standard e-equation is recovered. The Bardina model [7, 9] takes the form 


dv{ dvj 
dxj dxi t 


(10) 


e = c:, - c'J— 


K 


K 


(14) 


in a homogeneous turbulence, where 

I 

C* cl = 1.50 - 0.015^ 5 (15) 

CU = 1.83 + 0.15y (2^%) ’ (16) 

axe functions of the mean rotation tensor. The models (11), (12), and (14)-(16) are written 
for an inertial frame. In a non-inertial frame, the rotation tensor JDij must be replaced by 
the absolute mean rotation tensor 


Wij = u>ij + e mji Q m (17) 

where f l m is the angular velocity of the non-inertial frame relative to an inertial framing (see 
Speziale [10]). 
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Raj [8] developed a model based on an analysis of the exact transport equation of the true 
dissipation. In an inhomogeneous turbulence, this exact transport equation contains terms 
that depend explicitly on the rotation rate of the reference frame. Based on an empirical 
argument, Raj [8] developed a simple correction term to account for the effect of rotations 
which was guided by his turbomachinery research. For homogeneous turbulence, the model 
proposed by Raj [8] takes the form 


* — ~ + C t3 Qe ~ l) (18) 

in a rotating turbulent flow, where fi is the angular velocity of the reference frame - which, 
for simplicity, is aligned along the z direction - and t zx is the component of Reynolds stress 
tensor along the axis of rotation. Here, C c3 is a constant of 0(1). In the limited practical 
applications where this model has been used, C e3 has assumed different values ranging from 
1 to 5. For simplicity, we will set C e3 = 5 (the largest of these values); (7 el and C e2 assume 
the traditional values of 1.44 and 1.92, respectively. It must be said at the outset that (18) 
does not constitute a general model and has some internal inconsistencies in the way it was 
derived [11]. Nonetheless, the model does have some interesting features that are worth 
testing in the problems to be considered in this study (e.g., the model directly accounts for 
the effect of the rotationally induced anisotropy in the Reynolds stress tensor, along the axis 
of rotation, on the dissipation). 

In order to test these models in the flows to be considered, a Reynolds stress model is 
needed. A simple second-order closure model is chosen - namely, the shortened form of the 
Launder, Reece, and Rodi model which is now referred to as the “Basic Model” by Launder 
and co-workers [12]. For a homogeneous turbulence in a rotating frame, this Basic Model of 
Launder and co-workers takes the form [13]: 


T H 


= (C 2 ~ 1) T* 


&Uj 

dx k 


+ T jk 


dvA 

dx k ) 


-j-(Cj 2 )(Ti k e m kj^m T j k e mk i$l m ) 


(19) 


where C\ and C 2 are dimensionless constants that assume the values of 1.8 and 0.6, respec- 
tively. This model is derived by making two major assumptions - namely that the dissipation 
rate tensor is isotropic, i.e., 

= | ehj (20) 

and the pressure strain correlation is of the general form 

n a = eAi(b) + KMyui b)|| (21) 

where Aij and At tjk i are linear functions of the anisotropy tensor b. 
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The first problem to be considered as a test case is isotropic turbulence in a rotating 
frame. Here, an initially decaying isotropic turbulence wherein 

2 . 

T{j = ~K(}8xj } £ ~ £o> (22) 

at time t = 0, is subjected to a solid body rotation with constant angular velocity 

= (23) 


and vanishing mean velocity gradients 


dvj 

dxj 


= 0 . 


(24) 


The substitution of (23) - (24) into (19) - along with the use of the initial condition of 
isotropy (22) - then yields the equation 

K = -e (25) 


where Tij = | KSij. This is the same form as its non-rotating counterpart. Hence, the 
Launder, Reece, and Rodi model predicts that an initially isotropic turbulence subjected to 
a rotation decays isotropically - a result that is consistent with the Navier-Stokes equations 
[14, 15]. For rotating isotropic turbulence, the modeled e- transport equation takes the 
standard form 

i = -C.4 (26) 

(where (7*2 = 1.92) for the Pope and the Raj models; consequently, these models do not 
distinguish the difference between isotropic turbulence in a rotating frame and in an inertial 
frame. The Bardina model, however, does have a non-zero rotational correction; for rotating 
isotropic turbulence (14) takes the form 

e = -C e2 ^ - C e3 Ue (27) 

where (7*2 = 1.83 and (7*3 = 0.15. Likewise, so does the Hanjalic and Launder model for 
which we have the following e-transport equation in rotating isotropic turbulence: 

t = -<7*2^ - 4C e3 Kn 2 (28) 


where (7*2 = 1.92 and C* 3 = 0.27. 

The second problem to be considered is homogeneous shear flow in a rotating frame. The 
mean velocity gradients here take the form 


p- 0 S 0 \ 

-^= 0 0 0 (29) 

OXj 0 0 0 ) 
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which correspond to a plane shear. They are applied in a steadily rotating frame where the 
axis of rotation is normal to the plane of the shear, i.e., 

= (0,0, Q). 

The substitution of (28) - (29) into (19) yields the Reynolds stress transport equation 

Tij — TikAjk + TjkAik — 1.8— ^Tij — —K 6^ + — — (30) 

where 

■ 0 -fs + ffi 0 1 

Aij = -ffi 0 0 (31) 

.000 

and V = — T 12 S is the turbulence production. The e-transport equation for the Pope model 
is given by 

£ € ^ 

£ — ~C e \—Ti2S — C e 2 (32) 

in rotating shear flow. This is the same as the standard model due to the fact that the mean 
vortex stretching term in ( 11 ) vanishes for any plane homogeneous turbulence. The Bardina 
model for rotating shear flow is obtained by replacing C e \ and C e2 in (32) with 

C* tX = 1.50 - 0.015|iS - fi| K/e ( 33 ) 

c; 2 = 1.83 + 0.15|§5 - n\K/e. ( 34 ) 

Both the Hanjalic' and Launder model and the Raj model take the form 

5 6 ^ 

^ = —Oei—TuS — C e 2~^ + R (35) 

in rotating shear flow where C e i = 1.44 and C e2 = 1.92. For the Hanjalic and Launder model 

R = -C e3 K(S - 2 fl ) 2 (36) 

where C e 3 = 0.27, whereas for the Raj model 

R = c,3fi£ iw - 0 (37) 

where C e3 = 5. 

In rotating shear flow, the transport equations for T tJ and e are solved subject to the 
initial conditions of isotropy ( 22 ) at time t — 0 . 

3. Discussion of the Results 

For the two homogeneous turbulence problems to be considered, the models give rise 
to a coupled set of first-order nonlinear ordinary differential equations which are solved by 
a fourth-order accurate Runge-Kutta numerical integration scheme. First, the results for 
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isotropic turbulence in a rotating frame will be presented. For this problem, the solutions 
depend on Cl,K 0 , and e 0 only through the dimensionless parameter CIKq/cq. In Figures 
l( a ) - (b) the time evolution of the dimensionless turbulent kinetic energy (K* = K/Kq 
and r = eot/ K q is the dimensionless time) obtained for the Bardina model is shown for two 
moderate rotation rates: ClK 0 /e 0 = 0.123 and ClK 0 /e 0 = 0.469. These two rotation rates 
correspond to the test cases considered in the physical experiments of Wigeland and Nagib 

[15] . It is clear that the agreement between the model predictions and the experiments is good 
(however, it should be noted that the difference between these results and the predictions of 
the standard model are small). This is not surprising since the Bardina model was calibrated 
based on the same experimental data [7]. While the Bardina model performs well for weak 
rotation rates, it has substantial problems for rapidly rotating isotropic turbulence. In Figure 
1(c), the time evolution of the turbulent kinetic energy predicted by the Bardina model for 
ClK 0 /e 0 = 69.5 is compared with the direct simulations of Speziale, Mansour, and Rogallo 
[14]. It is clear that the model drastically underpredicts the decay rate of the turbulent 
kinetic energy. This arises from a deficiency in the model; it is a simple matter to show that, 
for QKo/cq > 1.11, the Bardina model predicts that the turbulence Reynolds number Re ( 
increases during the decay - a result that is unphysical. 

The time evolution of the turbulent kinetic energy predicted by the Hanjalic and Launder 
model is shown in Figure 2(a) along with the predictions of the standard model and the 
experimental data of Wigeland and Nagib on rotating isotropic turbulence corresponding 
to ClKo/e 0 = 0.123. While the model does not do too bad a job for this case, it yields 
worse results than the standard model. It then breaks down entirely for appreciable values 
of ClK 0 /e 0 . In Figure 2(b), the predictions of the Hanjalic and Launder model are shown 
for ClK 0 /e 0 = 0.469. The model yields extremely bad results for this case and becomes 
unrealizable when r = 4.2 (i.e., the dissipation rate becomes negative causing us to terminate 
the calculations). It is the source term added by Hanjalic and Launder to the right-hand-side 
of the dissipation rate equation that causes the model to become unrealizable for discernible 
values of CIKq/cq. We will not show the results for the Pope model and the Raj model since 
they are identical to the predictions of the standard model. Hence, like the standard model, 
they are unable to predict the reduction in the decay rate of the turbulent kinetic energy 
that arises from a system rotation of isotropic turbulence [14, 15]. 

Now, we will discuss the results obtained for rotating shear flow. In Figures 3(a) - (c), the 
time evolution of the turbulent kinetic energy predicted by the Bardina model is compared 
with the standard model and the large-eddy simulations of Bardina, Ferziger, and Reynolds 

[16] for three rotation rates: Cl/S = 0, 0.25, and 0.5. (Again, K* = K/K 0 whereas t * = St). 
From Figures 3(b) - 3(c) it is clear that the Bardina model has little effect for appreciable 
rotation rates; like the standard model, the Bardina model drastically underpredicts the 
growth rate of the turbulent kinetic energy for Cl/S = 0.25 and Cl/S = 0.5. However, what 
is more serious is that the Bardina modification degrades considerably the results for pure 
shear flow (Clj S — 0) as shown in Figure 3(a). The growth rate of the turbulent kinetic 
energy predicted by the Bardina model is too high. This results from the overprediction of 
the equilibrium value of the shear anisotropy (612)00 = -0.191; the standard model yields a 
value of (613)00 = —0.18 which is closer to the experimental value of (6 i 2 )oo = —0.16. 

In Figure 4, the predictions of the Hanjalic and Launder model are compared with the 
standard model and the large-eddy simulation of Bardina, Ferziger, and Reynolds [16] for 
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homogeneous shear flow (Cl/ S = 0). It is clear that the Hanjalic and Launder model predicts 
much too strong a growth rate for the turbulent kinetic energy in homogeneous shear flow. 
This results from an underprediction of the dissipation rate. In fact, the calculation is 
terminated a short time after St = 2 due to the model becoming unrealizable through 
negative dissipation rates. Since the Hanjalic and Launder model yields anomalous results 
for pure shear flow (the most basic case) we will not bother to show the results for non-zero 
rotation rates. 

In Figures 5(a) - (b), the results predicted by the Raj model for rotating shear flow are 
compared with those of the standard model as well as with the large-eddy simulations of 
Bardina, Ferziger, and Reynolds [16]. We do not show the results of the Raj model for pure 
shear flow (Cl/S — 0) since they are identical to those of the standard model. From Figure 
5(a), it is clear that the Raj model yields improved predictions over the standard model for 
the Cl/S — 0.25 case. However, there are problems with the Raj model for rotation rates 
that are appreciably larger than Cl/S = 0.25 as can be seen in Figure 5(b). For the case 
where Cl/S = 0.5, the Raj model yields unrealizable results (i.e., the dissipation rate becomes 
negative for values of St appreciably larger than 7). 

It is clear from these results that simple corrections to the modeled dissipation rate 
equation obtained by adding a source term that depends on rotational strains can give rise to 
realizability problems (see Lumley [1]). Since the standard modeled dissipation rate equation 
is realizable with respect to K and e (i.e., as shown by Speziale [17], the standard model will 
always yield positive kinetic energies and dissipation rates in homogeneous turbulent flows), 
it is clear that these kind of ad hoc corrections are counterproductive. In order to properly 
describe the turbulent flows discussed in this section, directional and two-point information 
is needed in the modeling of the dissipation rate. Speziale and Gatski [18] have been recently 
working on the development of a modeled transport equation for the tensor dissipation rate 


£ij 


= 2 v 


dui duj 
dx k dx k 


This model for the tensor dissipation takes the form 




Sik 


(fj + 2emkjUm 



dvi 

dx k 


+ 2 e m ki$l 



+ -C , e i— — C e i—£ij + C c $eSij 


(38) 


— C ei (£ikWjk + EjfcW'ifc) 

in a rotating frame where C e \ = 1.46, C e 2 — 1-83, C e 3 = 1.50, and C c 4 = 1.0. It was 
developed for use in conjunction with the SSG model for the pressure-strain correlation 
derived recently by Speziale, Sarkar, and Gatski [19]. In Figures 6(a) - (c), the results 
predicted by the tensor dissipation model (used in conjunction with the SSG model) are 
compared with the predictions of the Launder, Reece, and Rodi model and the large-eddy 
simulations [16]. It is clear that this new tensor dissipation model does a much better 
job in capturing the trends of the large-eddy simulations. Most notably, the substantially 
stronger growth rate of the tl/S = 0.25 case is, for the most part, captured; furthermore, 
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the model does not restabilize for the Cl/ S = 0.5 case - a result that is consistent with linear 
stability theory as well as with the large-eddy simulations. These improvements are not that 
surprising in light of the recent findings of Durbin and Speziale [20] which suggest that the 
dissipation rate tensor can be anisotropic in equilibrium homogeneous shear flows at high 
turbulence Reynolds numbers. 

4. Conclusions 

Several modifications to the modeled dissipation rate transport equation that have been 
proposed to account for rotational strains are compared critically for two test flows. The 
following general conclusions can be drawn: 

(1) The Pope [5] model yields the same predictions as the standard model for rotating 
isotropic turbulence and for rotating homogeneous shear flow. These predictions are deficient 
in that they fail to properly account for the reduction in the turbulence decay rate that results 
from a system rotation of isotropic turbulence and they fail to account accurately for the 
changes in the growth rate of the turbulent kinetic energy that result from a system rotation 
of homogeneous shear flow which can either stabilize or destabilize the flow. 

(2) The Hanjalic and Launder [6] model has major problems with the violation of real- 
izability. It yields negative dissipation rates in rotating isotropic turbulence for appreciable 
values of CIKq/cq. Furthermore, this modification degrades significantly the predictions for 
homogeneous shear flow (the predicted growth rate of the turbulent kinetic energy is much 
too large). 

(3) The Bardina model [7] is only able to accurately predict the reduction in the decay rate 
of the turbulent kinetic energy in rotating isotropic turbulence for weak to moderate rotation 
rates where this effect is small. For CIKq/eq 1, this model predicts too strong a reduction 
in the turbulence decay rate. Furthermore, this model does not yield any improvements for 
rotating homogeneous shear flow and mildly degrades the predictions for pure shear flow. 

(4) The Raj [8] model yields the same deficient predictions as the standard model in 
rotating isotropic turbulence. In rotating homogeneous shear flow, this model does yield 
some improvements for moderate rotation rates for which Cl/S < 0.3. However, for stronger 
rotation rates the model becomes unrealizable and yields negative values for the turbulent 
dissipation rate. 

In order to properly describe rotating turbulent flows, some directional and two-point 
information needs to be incorporated into models for the dissipation rate rather than the ad 
hoc modifications discussed above. The directional information can be incorporated through 
the use of the tensor dissipation. In this regard, the preliminary calculations presented in 
Section 3 for rotating shear flow based on a modeled tensor dissipation rate equation appear 
to be promising. However, the reduction in the decay rate of the turbulent kinetic energy 
in rotating isotropic turbulence is a two-point phenomena; the inertial waves generated by a 
system rotation scramble the transfer term in such a way that the phase coherence needed 
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to cascade energy from large to small scales is disturbed. Some limited two-point informa- 
tion needs to be included based on an appropriate integral length scale which responds to 
rotational strains. This issue is currently under investigation and will be the subject of a 
future paper. 
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Rotating Isotropic Decay 


(a) n Kq/s 0 = 0.123 


Figure 1. The decay of turbulent kinetic energy in rotating isotropic turbulence: 

Bardina model; standard model; o experimental data [15]; • simulations [14]. (a) 

SlKo/eo = 0.123, (b) QX 0 /e 0 = 0.469, (c) Q K 0 /e 0 = 69.5. 
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Figure 1. The decay of turbulent kinetic energy in rotating isotropic turbulence: 

Bardina model; — standard model; o experimental data [15]; • simulations [14]. (a) 

nKo/eo = 0.123, (b) UK 0 /e 0 = 0.469, (c) SlK 0 /eo = 69.5. 
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Figure 1. The decay of turbulent kinetic energy in rotating isotropic turbulence: 

Bardina model; standard model; o experimental data [15]; • simulations [14]. (a) 

QKo/eo = 0.123, (b) QK 0 /e 0 = 0.469, (c) SlKo/e o = 69.5. 






Rotating Isotropic Decay 


(a) Q K n /t 


Figure 2. The decay of turbulent kinetic energy in rotating isotropic turbulence: — 

Hanjalic and Launder model; standard model; o experimental data [15]. (a) QKo/en 

0.123, (b) nK 0 /eo = 0.469. 
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Figure 2. The decay of turbulent kinetic energy in rotating isotropic turbulence: — 

Hanjalic and Launder model; standard model; o experimental data [15]. (a) Q,K 0 /e 0 

0.123, (b) ClKo/eo = 0.469. 
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Figure 3. Time evolution of the turbulent kinetic energy in rotating shear flow for e 0 / SK 0 — 

0.296: standard model; Bardina model; o large eddy simulations [16]. (a) Cl/S = 0, 

(b) Cl/S = 0.25, (c) Cl/S = 0.5. 





Figure 3. Time evolution of the turbulent kinetic energy in rotating shear flow for e 0 /SKo = 

0.296: standard model; Bardina model; o large eddy simulations [16]. (a) f2/5 = 0, 

(b) Q/S = 0.25, (c) fi/S = 0.5. 
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Figure 4. Time evolution of the turbulent kinetic energy in homogeneous shear flow for 

eo/SKo = 0.296: standard model; Hanjalic and Launder model; o large eddy 

simulation [16]. 
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Figure 5. Time evolution of the turbulent kinetic energy in rotating shear flow for e 0 /SK 0 = 

0.296: standard model; Raj model; o large eddy simulations [161. fa) Cl/S = 0 25 

(b) Cl/S = 0.5. 




Figure 5. Time evolution of the turbulent kinetic energy in rotating shear flow for eo/SKo = 

0.296: standard model; Raj model; o large eddy simulations [16]. (a) Q/S = 0.25 

(b) Cl/S = 0.5. 
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Figure 6. Time evolution of the turbulent kinetic energy in rotating shear flow for e 0 /SK 0 = 

0.296: standard model; tensor dissipation model [181; o large eddy simulations 

[16]. (a) Cl/S = 0, (b) Cl/S = 0.25, (c) Cl/S = 0.5. 
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Figure 6. Time evolution of the turbulent kinetic energy in rotating shear flow for e 0 /SK 0 = 

0.296: standard model; tensor dissipation model [18]; o large eddy simulations 

[16]. (a) Q/S = 0, (b) a/S = 0.25, (c) Q/S = 0.5. 
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Figure 6. Time evolution of the turbulent kinetic energy in rotating shear flow for e 0 /SK 0 = 
0 296: standard model; tensor dissipation model [18]; o large eddy simulations 

[16]. (a) a/s = o, (b) a/s = o.25 ( (c) a/s = 0.5. 
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